Velocities as a probe of dark sector interactions 



o : 
o . 

> ' 

\o . 

u: 
43 : 

Oh- 
i < 

o ■ 



> 

(N 

o 

O 



13 



Kazuya Koyama, Roy Maartens, Yong-Seon Song 
Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth P01 3FX, UK 

(Dated: October 26, 2009) 

Dark energy in General Relativity is typically non-interacting with other matter. However, it is 
possible that the dark energy interacts with the dark matter, and in this case, the dark matter can 
violate the universality of free fall (the weak equivalence principle) . We show that some forms of the 
dark sector interaction do not violate weak equivalence. For those interactions that do violate weak 
equivalence, there are no available laboratory experiments to probe this violation for dark matter. 
But cosmology provides a test for violations of the equivalence principle between dark matter and 
baryons - via a test for consistency of the observed galaxy velocities with the Euler equation. 



I. INTRODUCTION 

Dark matter is currently only detected via its gravi- 
tational effects, and there is an unavoidable degeneracy 
between dark matter and dark energy within General 
Relativity. There could be a hidden non-gravitational 
coupling between dark matter and dark energy, and thus 
it is interesting to develop ways of testing for such an 
interaction (see [l|, 0, 9 for earlier attempts). 

One signal of a dark sector interaction could be a viola- 
tion of the weak equivalence principle (universality of free 
fall) by dark matter, under the non-gravitational drag 
due to coupled dark energy. Since Galileo shattered the 
myth that heavier objects fall faster, the universality of 
free fall has been established as a fundamental principle 
of gravity. Laboratory tests have been made to show the 
independence of the acceleration of objects from their 
masses and chemical composition. However, these tests 
apply to baryonic matter, and no direct probe of dark 
matter is available. 

If the interacting dark sector couples non- 
gravitationally to baryonic matter, then existing 
laboratory tests provide constraints on the dark sector 
interaction Q. Here we assume that there is zero (or 
negligible) non-gravitational coupling between the dark 
sector and standard-model fields. A difference in the 
acceleration between dark matter and baryons could 
show up in the stellar distribution in tidal trails of satel- 
lite galaxies [f|. This same difference should also show 
up as an inconsistency when interpreting the relation 
between galaxy peculiar velocities and overdensities, as 
we explain below. 

We assume that gravity on all scales is described by 
General Relativity. Thus there is no gravitational mech- 
anism to violate the weak equivalence principle. Note 
that this is also true of scalar-tensor theories, since the 
gravitational scalar degree of freedom couples equally to 
all types of matter. Indeed, most metric theories of mod- 
ified gravity also respect the weak equivalence principle 
(see, e.g., |g|). Various tests have been developed to dis- 
criminate between metric theories of modified gravity, 
and non-interacting dark energy models in General Rel- 
ativity (see, e.g., 0]). But these tests do not in general 
apply to dark energy that interacts with dark matter, 



since a dark sector interaction can introduce new degen- 
eracies [||. We confine ourselves to the question of how 
galaxy peculiar velocities can be used to detect dark sec- 
tor interactions within General Relativity. 



II. INTERACTING DARK ENERGY 

We briefly review the necessary background on per- 
turbations of interacting dark energy models in General 
Relativity. (For recent work with further references, see, 
e.g., El.) 

A general dark sector coupling may be described in the 
background by the energy balance equations of cold dark 
matter (c) and dark energy (x), 



p' c = -3TLp c + aQ c , 



P.r 



-3H(l + w x )p x + aQ x , Q x = -Q c 



(1) 
(2) 



where 



time, with ds 2 



Px/Px, "H = dlna/dr and r is conformal 



-dr + dx ). Here Q Cl Q x are the 



rates of energy density transfer to dark matter and en- 
ergy respectively. In order to avoid stringent "fifth-force" 
constraints, we assume that baryons (6), photons (7) and 
neutrinos (y) are not coupled to dark energy and are sep- 
arately conserved. 

In the Newtonian gauge the perturbed metric is given 
by 



ds 1 



(1 + 2^)dr 2 + (1 - 2V)dx : 



(3) 



where we have neglected anisotropic stress since we are 
interested in the late universe. The total (energy-frame) 
four- velocity is 



(4) 



= a' 1 [l - #,<9V 
where the velocity potential v is defined by 
( P + P)v = J2(pa+Pa)va, 
and A = c, x, b, 7, v. The A-fluid four-velocity is 
u'X = a- 1 (l-y,d i v A K 



(5) 



(6) 
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The covariant form of energy-momentum transfer is 

V„Tf = Q\ , (7) 

where Q A = for A = 6, 7, v in the late universe, while 
QJf = — Q% ^ 0. The energy-momentum transfer four- 
vector can be split relative to the total four-velocity as 

Qa = Qau» + Qa = Qa + SQa , UpFZ = , (8) 

where Qa is the energy density transfer rate and F A is the 
momentum density transfer rate, relative to . Then it 
follows that F A = a -1 (0, d % Ja), where f A is a momentum 
transfer potential, and 



Qo = -a[Q A (l + y) + 5Q A ] 
Qf = adi(f A + Q A v). 



(9) 
(10) 



In the background, the energy-momentum transfer four- 
vectors have the form Q^ = a^ 1 (Q c ,0) = — Q£ , so that 
there is no momentum transfer. 

The evolution equations for the dimensionless density 
perturbation 5 A = 5 pa/ Pa and for the velocity pertur- 
bation are: 

5' A + m(c 2 sA - w A )S A - (1 + w A )k 2 v A 
- 3H[3H{1 + w A ){c 2 sA - w a) + w' a ]va 



-3(l + w A )V = — SQ A 

PA 

+ ?QA[*-6 A -3H(c 2 sA -w A )v A ] 

PA 

y A +n(i-3c 2 sA ) VA 



(ii) 



(1 + w A )p A 



(1 + wa) 
[Q A [v-{l + c 2 sA )v A ] +/a},(12) 



where w c = = c 2 sc and c 2 x = 1. 

For our purposes, we are interested in the behaviour 
of dark matter in the Newtonian regime on sub-Hubble 
scales. In this case, the perturbed continuity and Eulcr 
equations reduce to 



5' c -k 2 v c = —{5Q C -Q C 5 C ), 



(13) 



v' c + Hv c + * = —[Q c {v-v c )+f c ], (14) 

Pc 

If the right-hand side of the continuity equation (|13|) is 
nonzero, then the interaction will lead to a bias in the 
linear regime between dark matter and baryons [l2| , since 
the baryon overdensities obey 



(15) 



S' b — k 2 vi, = . 



If the right-hand side of the Euler equation (jl4|) 
is nonzero, then the dark matter no longer follows 
geodesies and breaks the weak equivalence principle, un- 
like baryons, for which 



In the Newtonian regime, the Poisson equation be- 
comes 



k 2 ^) = -AnGci 2 (p c S c + p b 5 b ) 



(17) 



Here we neglect dark energy clustering, assuming that 
the sound velocity of dark energy perturbations is c sx = 
1 . Dark energy perturbations can be important on large 
scales depending on the strength of interactions but they 
are not important on sub-horizon scales as long as the 
sound velocity of dark energy perturbations is positive 
- since in that case, the gradient term in the evolution 
equation for 5 X [see Eq. (fT2"|) ] always dominates over the 
interaction terms. The evolution equation for 5 C is then 
given by 



S" c + H5' c - 4irGa 2 Pc S c - H-(5Q C ~ Q C 5 C 

Pc 



— {5Qc ~ QcS c 

Pc 



--k 2 

Pc 



Qc{v - V c ) + f c 



= 0. 
(18) 



III. DIFFERENT TYPES OF INTERACTION 

There is no fundamental theory that determines the 
form of the interaction, i.e., of Q^, so we are forced to 
use phenomenological models. Here we consider three 
types of interaction, each illustrated with a particular 
form: interactions that do not change the continuity or 
Euler equations; interactions that change only the Eulcr 
equation; interactions that change only the continuity 
equation. The general case, where both equations are 
modified, can be thought of as a linear superposition of 
the last two cases. 



A. Continuity and Euler equations unchanged 

A general class of interactions may be defined by re- 
quiring that there is no momentum exchange in dark mat- 
ter rest frame, 



= Qc< , 



(19) 



where Q c remains to be specified. For this class, we find 
from Eqs. © and (fTU|) that, for any Q c , we have f c — 
Qc(v c — v). Thus Eq. (fl4|) becomes 



Hv c + * = . 



(20) 



This is the same Euler equation as the non-interacting 
case, so that the dark matter velocity is not directly 
affected by the interaction and there is no violation of 
weak equivalence. The dark matter continues to follow 
geodesies, and feels no direct drag force from the dark 
energy. 

An example in the form of Eq. (19]) is [|, 0, HI 



v' b + Hv b + * = . 



(16) 



-rp c < 



(21) 
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where T is a constant interaction rate. In this case Q c 
—Tp c (l + S c ) and Eq. (fTB")) becomes 



6' - fc V = . 



(22) 



The continuity equation is therefore the same as in the 
non-interacting case. 

Thus for this form of interaction, there is no violation 
of the weak equivalence principle by dark matter, and no 
bias is induced by the interaction. In fact, in the Newto- 
nian regime, the only signal of the dark sector interaction 
in structure formation to linear order is via the modifica- 
tion of the background expansion history. The evolution 
equation l|18p for S c becomes 



% + H8' c - A7rGa 2 (p c 5 c + Pb S b ) = 0, 



(23) 



which is the same as in the uncoupled case. Thus the 
only imprint of the dark sector interaction on 8 C is via 
the different background evolution of Ti and p c . 



B. Continuity equation modified 

If we keep Eq. (Q3JJ) but generalize Eq. ([ST]) to [l5j 

Q£ = -(r c p c + r xPx )<, (24) 

then SQ C — Q C S C — T x p x (6 c — S x )- Since dark energy does 
not cluster on sub-Hubble scales, we can neglect the S x 
term, and we have 



5' — k 2 v r = aT. 



Px 



(25) 



For this interaction, the dark matter continues to follow 
geodesies by virtue of Eq. (|2"0|) , but the continuity equa- 
tion (|25l is modified. As a consequence, there will be a 
bias induced by the interaction. 

The evolution equation (fTg)) for 5 C becomes 



5" 



U-aT x ^- ) S' c = 4nGa 2 p b 6 b 



AirGa 2 p c + 2aHT x — 

Pc 



aT 



■(£)' 



(26) 



(This generalizes [15j , where only the case r c = is con- 
sidered.) 

The modification of the standard evolution for S c oc- 
curs in 3 ways: firstly via the modified expansion history 
in the background Ti and p c ; secondly by the modified 
Hubble friction term Ti —> Ti[l—aT x p x /Tip c ]', and thirdly 
by the modified effective gravitational coupling for dark 
matter - dark matter particle interactions, 



G e g = G 



1 



Up x 

2irGap2 4irGap c V p, 



+ 



x (Px\' 

'aOr \ On' 



(27) 



C. Euler equation modified 

A second general class of interactions has no momen- 
tum exchange in the dark energy frame, 

Q£ = Q c < . (28) 
It follows that f c = Q c {v x — v), and hence 

v' c + Tiv c + V=— Qc(vx-v c ). (29) 

Pc 

In this case, there is an explicit deviation of the dark 
matter velocity relative to the non-interacting case. The 
dark matter no longer follows geodesies in general. Note 
that, even though dark energy does not cluster on sub- 
Hubble scales, we cannot in general neglect the dark en- 
ergy velocity v x relative to the dark matter velocity v c in 
Eq. (EU). 

An example of the form of Eq. ([28]) is 

(30) 



Q£ = -ap c VV, 



where ip is the scalar field that describes dark energy and 
a is a coupling constant. Note that V M </? is parallel to 
the dark energy four- velocity u^: 



a V <p' 



(3D 



In this case, Q c — a 1 a{p c Lp' + 8p c ip' + p c 8(fi' — Pcf'^)- 
The perturbed Klein-Gordon equation is [lTf 

Sip" + 2W<V + (fc 2 + a 2 V w )5 V 

= 2<p'(^' + TW) + 2ip"V -aa 2 p c S c , (32) 

where V((p) is the quintessence potential. In the Newto- 
nian regime, the last term on the right dominates over 
the other terms, while the k 2 term dominates on the left, 
leading to 

k 2 5ip = —aa 2 p c S c . (33) 
It follows from Elqs. J29]), d30j) and fl3I) that 

v' c + Tiv c + y = -atp'(v c + ^y), (34) 

confirming the violation of weak equivalence. For the 
perturbed continuity equation (|13p . the right-hand side 
becomes —aSip'. By Eq. (f3"3")> . this term is suppressed 
by fc~ 2 relative to the S' c term on the left-hand side, and 
therefore to a good approximation we have 



5' - fc V = . 



(35) 



Using (f3"3"| and (|3"5|) , the evolution equation (fT5|) for S c 

becomes 

8" c +(Ti + a<p'% = 4irGa 2 p b 6 b 
+ 4nGa 2 (l + £^)p c S c . (36) 



4 



As in the case of Eq. (|26|) . the modification of the 
standard evolution for S c occurs in 3 ways [Hj]: firstly 
via the modified expansion history in the background H 
and p c ; secondly by the modified Hubble friction term 
H — ► H[l + atp'/Ti]; and thirdly by the modified effec- 
tive gravitational coupling for dark matter - dark matter 
particle interactions, 

Goff = G(l + ^). (37) 

These effects are incorporated in the modified iV- body 
simulations for this form of interacting dark energy [18j |. 



IV. TESTING FOR DARK SECTOR 
INTERACTIONS 

In this section, we discuss several possible ways to use 
observations to constrain the dark sector interactions dis- 
cussed in the previous section. 
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FIG. 1: The ratio between the true matter density obtained 
from structure formation and the density estimated from ge- 
ometrical tests assuming the non-interacting adiabatic be- 
haviour p m oc a -3 . 



A. Continuity and Euler equations unchanged 

We first consider the case where there is no modifi- 
cation to the dynamics of perturbations in the Newto- 
nian regime. The difference comes purely from the mod- 
ified background evolution. If dark matter interacts with 
dark energy, the dark matter density no longer decays 
like a -3 . This affects the distance measures in the Uni- 
verse and thus changes the measurements of CMB, SNe 
and Baryon Acoustic Oscillations. By combining these 
observations, we can measure today's matter density and 
then determine the matter energy density at the last scat- 
tering surface. However, the distance is determined by 
integrating over the expansion history and we cannot di- 
rectly check the deviation at each redshift from the stan- 
dard behaviour, p c oc a~ 3 . 

There is an independent way to measure the dark mat- 
ter density using structure formation. From the Poisson 
equation, the dark matter density can be written as 



w m (o) = Vt m {a)h 2 = - 



2*0, a) / kh 



35 c (k,a) \aHo 



(38) 



where we neglected the baryon contribution for simplic- 
ity (we are only illustrating the principle, rather than 
making quantitative predictions). One way to measure 
5 C is to reconstruct 6 C from peculiar velocities using the 
continuity equation Eq. (|22p because in this case there is 
no modification to the continuity equation and no viola- 
tion of weak equivalence principle. On the other hand, 
weak lensing measures directly * without bias. Thus we 
can use Eq. (|38[) to predict the background evolution of 
matter density from structure formation. 

In Fig.[TJ we plot w m /w£f, where iv m is the true mat- 
ter density measured by weak lensing, and w^f is derived 



from the background measurement of u> m at the last scat- 
tering surface, assuming p m oc a~ 3 . At late times when 
the interaction becomes important, this ratio deviates 
from I due to the non-adiabatic decay of the dark mat- 
ter density. In this way, we can check the modification 
to the behaviour of the matter density at each redshift, 
using tomographic measurements of W from weak lensing. 



B. Test of the continuity equation 

In the case where the interaction changes only the con- 
tinuity equation, there is no difference between the pe- 
culiar velocities of baryons and dark matter. We assume 
that galaxies can be treated as test particles that are 
made of baryons and whose peculiar velocities, v g , are de- 
termined by baryon peculiar velocities. Although there 
is an indication that this assumption is valid [19j], this 
should be tested by N-body simulations carefully in the 
presence of interaction. We leave this for a future work. 

With this assumption, we can determine peculiar ve- 
locities of baryons, Vb, from peculiar velocities of galaxies, 
v g . The latter can be measured for example by redshift- 
space distortions (see [2(J Hl[ for recent work) . Then it 
is possible to determine dark matter peculiar velocities 
because v c = U&. 

On the other hand, density perturbations can be mea- 
sured from the galaxy distribution with a knowledge of 
bias. One possibility to measure bias is to use weak lens- 
ing. Weak lensing measures ^> without bias and S c can 
be derived from the Poisson equation (| 1 T|) . Note that in 
order to measure \& from 8 C , it is necessary to measure 
the true evolution of p Cl which is modified by interac- 
tions. However, we found that the modification to the 
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case. Thus, if we use the estimated peculiar velocity, the 
continuity equation is apparently broken 



„2„,cst 



(40) 



In this case, the continuity equation is not broken but 
Vc 7^ v b- Then we can apply the same analysis as in the 
previous section. We can measure 5 C from weak lensing. 
Then it is possible to prove the breakdown of the weak 
equivalence principle through the apparent breakdown of 
the continuity equation. Fig. [3] demonstrates this idea. 



V. CONCLUSIONS 



FIG. 2: The breakdown of the continuity equation by an in- 
teracting dark energy model. In this model, Vb = v c and 
S' c - k 2 v c = aT x {p x / p c )5 c . 

continuity equation has significant effects even in weak 
interactions cases where the effect of interactions on p c 
is negligible. Thus in the following we only consider the 
case where we can neglect the effect of interactions on the 
evolution of p c . Another possibility is to use the pecu- 
liar velocity measurements. In the case that we consider 
here, the Euler equation is not modified [see Eq. ([20]) ] 
and it is possible to reconstruct \& from v c . Then again 
using the Poisson equation, S c can be derived [22II23II24I]. 
In this way we can test whether the continuity equation 
is modified. 

Fig. [5] demonstrates the breakdown of the standard 
continuity equation by an interacting dark energy model. 
We used a model where T x ^ and r c = 0. 

C. Test of weak equivalence principle 

The weak equivalence principle is broken when the 
Euler equation for dark matter is modified. In this 
case, there is a difference between the peculiar veloci- 
ties of dark matter and baryons. With the assumption 
that galaxies trace baryon peculiar velocities, we measure 
baryon peculiar velocities say from red-shift distortions. 
Unlike the previous case, dark matter peculiar velocities 
are different. However, without knowing that there is 
an interaction between dark matter and dark energy, we 
estimate dark matter peculiar velocity as 

vT =v b =v g . (39) 

The estimated peculiar velocity is different from true pe- 
culiar velocity of dark matter: v c ^ v° st . 

If the continuity equation is not modified, as happens 
for the model of Eq. (|3H)l , then the true peculiar velocity 
satisfies the same continuity equation as the uncoupled 



An interaction between dark matter and dark energy 
could exist in various ways which are not detectable by 
any direct probe. We investigated how the Euler equa- 
tion and the continuity equation for dark matter could be 
modified by such an interaction, taking care to provide 
a covariant analysis of momentum transfer. Modifica- 
tion of the Euler equation indicates a deviation of the 
dark matter motion from geodesic, under the drag force 
of dark energy - and a consequent breaking of the weak 
equivalence principle for dark matter. Using three dif- 
ferent forms of interaction as examples, we considered 
interacting models in which: 

(A) neither the Euler nor continuity equations are modi- 
fied, so that the effect of the interaction in the Newtonian 
regime is purely via the different background evolution; 

(B) the Euler equation is unchanged but the continuity 
equation is modified (and consequently a new bias is in- 
troduced by the interaction); 

(C) the continuity equation is unchanged but the Euler 
equation is modified, leading to violation of weak equiv- 
alence. 

We discussed how in principle observations could be 
used to detect these different forms of interacting dark 
energy. In case (A), we used the fact that the continuity 
and Euler equations are unchanged to devise a test based 
on the non-adiabatic redshifting of the dark matter. This 
test uses independent measurements of the Newtonian 
potential and the density perturbation via the Poisson 
equation, to compute the true matter density and show 
that it deviates from the non-interacting case. 

In cases (B) and (C), the effects of the violation of the 
continuity equation or Euler equation are stronger than 
the non-standard redshifting of the background matter 
density. We showed how, given a knowledge of bias from 
weak lensing, tests could be constructed for the break- 
down of the continuity or the Euler equation. 

A further issue raised by our investigations is how to 
distinguish interacting dark energy from modified grav- 
ity. This is left for future work. 
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